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New insights on the interior of solar-like pulsators 
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Abstract 

The high performance photometric data obtained 
with space mission CoRoT offer the opportunity to 
efficiently constrain our models for the stellar inte- 
rior of solar-like pulsating stars. On the occasion of 
the analysis of the oscillations of solar-like pulsator 
HD 49385, a GO-type star in an advanced stage of evo- 
lution , we revisit the phenome n on of the avoided cross- 
ings. IChristensen-Dalsgaard I (|2003l ) proposed a sim- 
ple analogy to describe an avoided crossing between 
two modes. We here present an extension of this anal- 
ogy to the case of n modes, and show that it should 
lead, in certain cases, to a characteristic behavior of the 
eigenfrequencies, significantly different from the n = 2 
case. This type of behavior seems to be observed in 
HD 49385, from which we infer that the star should be 
in a Post Main Sequence phase. 

Keywords Stellar evolution - Stellar pulsation - Stel- 
lar interior - Advanced stages of stellar evolution 



1 Introduction 

As a star evolves, the eigenfrequencies of the g modes 
increase due to an increase of the buoyancy frequency. 
When the frequency of a non-radial g mode becomes 
close to that of a p mode of same degree £, the two 
modes undergo an avoided crossing, at the end of which 
they have exchanged natures. During the avoided cross- 
ing, both modes present a mixed character: they act as 
g modes in the deep interior, and as v modes below th e 
surface (see e.g. IScuflairdll974l . lAizenman et al.lll977l ). 
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The avoided crossings happen on very short time- 
scales compared to the evolution timescale. Besides, 
the existence of mixed modes and their frequencies if 
they exist, depend greatly on the profile of the buoy- 
ancy frequency throughout the star. The observation 
of such modes would therefore yield strong constraints 
on the age of the star, and on its inner structu re, as 
described in Dziembowski fc Pamvatnvkhl ( 1991 ). For 
solar-like pulsators, though models showed that mixed 
modes should be present in t he spectrum of certai n 
evolved objects, e.g. r) Bootis ( Di Mauro et al. 2003| ). 
they were never clearly identified. 

Avoided crossings are usually assumed to involve 
two modes only. It is the case, for exemple, in 



Christensen-Dalsgaard I (|2003l ) , where a simple analogy 
is proposed to describe the phenomenon. We present in 
Sect. an extension of this analogy to the case of an 
avoided crossing between n modes and apply our results 
to study the evolutionary status of solar-like pulsator 
HD 49385 in Sect. El 



2 Avoided crossing with n modes 

2.1 A simple analogy for n = 2 

The fact that modes with close eigenfrequencies have a 
mixed character is due to the coupling which exists be- 
tween the g-mode cavity and the p-modc cavity. If they 
were uncoupled, two modes of different nature could 
share the same eigenfrequency, without any impact 
on their eigenfunctions . Based on an original idea of 
yon Neuman fc Wiener ( 1929I ). Christensen-Dalsgaard 



( 20031) proposed a simple analogy to understand the 



process of the avoided crossings. He considered the two 
cavities of the star as a system of two coupled oscillators 
yi(t) and ?/2(0 with a time dependence, responding to 
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the following system of equations: 

d 2 yi(t) 



2.5 



dt 2 

d 2 y 2 (t) 
dt 2 



^ 2 (A) 2 j/ 2 + "1,22/1 



(1) 
(2) 



where 01,2 is the coupling term between the two oscil- 
lators. cJi(A) and u>2(A) are the eigenfrequencies of the 
uncoupled oscillators (0:1,2 = 0). They depend on a pa- 
rameter A, used to model the change of the dimensions 
of the cavities as the star evolves. We suppose that for 
a certain A = Aq, the uncoupled oscillators cross, i.e. 
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wi(A ) = W2(A ) = 



(3) 



We aim at determining the cigenfrcquency uj(X) of the 
whole system. The solution of Eq. [1] and [5] is under the 
form: 



Fig. 1 Variations of the eigenfrequencies uj± of the system 
with parameter A, for different values of the coupling term 
0:1,2- 
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By inserting Eq. 0] and O into Eq. [Tj and [2J we find 
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that the eigenfrequencies of the system are obtained by 














solving the eigenvalue problem AC = lo 2 C ', where 
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and C = [ci , c 2 ] . 
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We get the two following solutions for the system: 



1 



±2 [( w i- w 2) 2 +4oi, 2 ; 



1/2 



(7) 



We verify here that if the coupling term oi, 2 is very 
small compared to the difference between the eigenfre- 



Fig. 2 Variations of the coefficients c\ and c\ with 
parameter A, for different values of 0:1,2 (the lines correspond 
to the same values of ai,a as in Fig. [TJ. We note that 
the coefficients corresponding to the other eigenfrequency 
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quencies (01,2 <C 



4\), then the eigenfrequencies c^ +) (A) and 4 _) '(A) = 4 +) (^))- 



(A) are obtained by interchanging ci and C2 (c[ (A) = 



of the system are close to uj\ and UJ2- If, on the con- 



trary, 



a 1,2, then the eigenfrequencies can 



be approximated by ui± 



± oi 2- We see here that 



the two oscillators "avoid" the frequency ujq, and the 
larger 01,2, the further away of loq they are during the 
avoided crossing. 

To analyze these solutions in more details, we im- 
pose the variations of the eigenfrequencies with A: 
wi(A) 2 = 1 and W2(A) 2 = A for example. In this 
case, we can plot the eigenfrequencies ^(A) (see Fig. 
[lj. Normalizing the coefficients in such a way that 
ci(A) 2 + 02(A) 2 = 1, we obtain the variations of the 
coefficients with A, shown in Fig. [2] We clearly see 
that the modes exchange natures during the avoided 
crossing. 

The choice we made for the profile of u>i (A) and 0J2 (A) 
is not an innocent one. Indeed, if we represent the mode 



frequencies normalized by y GM/R 3 , the p-mode fre- 
quencies are almost constant during an avoided cross- 
ing, while the g-modc frequencies increase. 

2.2 Extension of the analogy to the case n > 2 

The previous analogy relies on the fact that two modes 
only are affected during an avoided crossing: this cor- 
responds to neglecting the coupling between the two 
considered modes, and the other modes in the spec- 
trum. Let us now push the analogy a step further, and 
consider the case where these other coupling terms play 
a role. We consider n oscillators, instead of two only. 
We take n — 1 of them to have constant eigenfrequencies 
with A, uif — i, i = 1, ...,n — 1 (to simulate p modes), 
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Fig. 3 Variations of the eigenfrequencies of (n— 1) p modes 
coupled to a g mode which undergoes avoided crossings with 
the p modes (here, n = 10) . The dashed lines correspond to 
a "weak coupling" (a — 0.05), and the full lines to a "strong 
coupling" (a = 0.35). 



and one of them to variate as u\ = A (to simulate a g 
mode). We introduce a coupling between the g mode 
and all the p modes, resulting in the following system 
of equations: 



dt 2 



i 2 y n -i{t) 
dt 2 

d 2 y n (t) 
dt 2 



-wi(A) 2 yi + ax^Vr, 



(8) 



-w„_i(A) 2 2/„_i + a n _i, n y n 

-uj n (X) 2 y n + ai jn yi + . . . + a„-i,„2M-i 



For simplicity, we assume here that the term of coupling 
with the g mode is equal to the same value a for all the 
different p modes. By writing the different oscillators 
Ui{t) = Ci exp(— iuit), the eigenfrequencies of the system 
are once again found by solving the eigenvalue problem 
AC = ui 2 C with 
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and C = [ci, ■ • • , c„]. 

The solution is plotted in Fig. [3] for two different 
values of a. 

Once again, the choice of the evolution of the u>i with 
A was made on purpose. The n—1 modes simulating the 
p modes were taken equidistant, as are the p modes at 



Fig. 4 The eigenfrequencies for a certain value of A are 
presented in an echelle diagram. In a, we observe the per- 
turbation of the ridge due to an avoided crossing in the case 
of a weak coupling (a — 0.05). In b, the perturbation in the 
case of a strong coupling (a — 0.35). In both figures, the 
dashed line shows the ridge the n — 1 first oscillators would 
follow without avoided crossing (or during an avoided cross- 
ing with no coupling between the modes). 



the first order of the asymptotic theory. For each value 
of A, wc can plot an echelle diagram, which is a com- 
monly used tool to characterize equidistances. In the 
case of uncoupled oscillators, we obtain a straight ridge 
(see Fig. H} . Fig. 0] shows the impact of the avoided 
crossings on the curvature of the ridge: in the case of a 
weak coupling between the modes, the approximation 
that only two modes are affected is legitimate since the 
rest of the ridge is the same as for a = 0. However, as 
the coupling term increases, this ceases to be true, and 
we can observe a significative change in the curvature 
of the ridge. 

2.3 Interest in stellar modelling 

Our analogy suggests that, provided the coupling be- 
tween the p-mode cavity and the g-mode cavity is 
strong enough, an avoided crossing does not affect only 
two modes of same degree £, but also a certain number 
of neighboring modes of degree £ which act mainly as 
p modes, but also have a small g-mode behavior in the 
center. This results in a distortion of the ridge of degree 
£ such as that shown in Fig. (3}d. The existence of such 
a distortion in the spectrum of an evolved star would 
make it much easier to detect the presence of a mixed 
mode and could give constraints on the inner structure 
of the star, as will be stressed in Sect. H) 

The strength of the coupling between the two cavities 
is inversely proportional to the size of the evanescent 
zone which separates them. The coupling will therefore 
be strongest for modes of low degree £, since their Lamb 
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frequency Si is smaller. We therefore expect the distor- 
tion of the ridge to be most significant for modes of low 
I. This is interesting since the modes of low degree are 
those we are most likely to observe in stars. We also 
remark that, since avoided crossings occur later in the 
evolution of the star for modes of lower £, we expect to 
observe this distortion only for very evolved objects. 



3 Application to HD 49385 

HD 49385 is a solar-like pulsator. of type GO and was 
observed with space mission C oRoT during 13 7 days . 



Through a detailed analysis, iDeheuvels et al. (|2009h 
found for the star T oS = 6095 ±50 K, logs = 4.00 ±0.05 
and a metallicity of [Fc/H]=0.09 ± 0.05 dex. The 
solar-like oscilla t ions of HD 49385 were analyzed in 
Deheuvels et al. ( 2009f ). An extensive modelling of the 
star is out of scope here, and will be presented in a 
future work. Our aim here is not to find an optimal 
model, but to show that the phenomenon we described 
in Sect. [5] seems to occur in HD 49385, and that it gives 
precious information regarding the evolutionary status 
of the star. 

3.1 Properties of the model 

Models were comput ed with the evolution code CE- 
SAM2k (jMorell 19971 ). and the mode frequencies were 
derived from the models using the L iege Oscillation 
Code (LOSC, see IScuflaire et alJl2008h . 

We used, for all our models, the OPAL 2001 equa- 
tion of state, and opacity tables as described in 



Lebreton et al.l (|2008l ). The nucle ar rates are compute d 



using the NACRE compilation f Angu lo et alJ 119991). 
Our m odels use the abundances of iGrevesse fc Noels 
(1993). The atmosphere is described by Eddington's 
grey temperature - optical depth law. The convective 
regions are treated usin g the Canuto, Goldm an, Mazz- 
itelli (CGM) formalism ( Canuto et al. 19961 ) calibrated 
on the Sun. 

For this preliminary modeling, we considered mod- 
els not including microscopic diffusion. An arbitrary 
amount of overshooting (extension of the convective 
core over a fraction a ov = 0.1 of the pressure height 
scale) was added. 

3.2 Evolutionary stage of HD 49385 

Fig. [5] shows that we can find models both in the Main 
Sequence (MS) and Post Main Sequence (PoMS), which 
fit the position of HD 49385 in the HR diagram. We are 
therefore facing a problem which is quite classical for 
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Fig. 5 Two evolutionary tracks fitting the position of 
HD 49385 in the HR diagram. The stellar parameters of 
HD 49385 are indicated within 1-cr error bars by the box. 
The full line corresponds to a MS model, and the dashed 
line, to a PoMS model. 



evolved object: a degeneracy of MS models and PoMS 
models which results in an uncertainty on the evolu- 
ti qnary stage of the st udied obj ect (see e.g. Procyon A 



Provost et alll2006l . 77 Boo in lDi Mauro et al.ll2003h . 



The use of a mean value of the large spacing does not 
cither provide enough constraint on the star to discrim- 
inate between the two stages of evolution. 

To investigate further on this problem, one must 
make use of the fact that MS models and PoMS mod- 
els fitting the position of HD 49385 in the HR diagram 
have a very different structure in their deep interior. In- 
deed, MS models are still burning hydrogen in the core, 
which is convective, whereas PoMS models are burning 
hydrogen in shell. Their core is no longer convective, 
and almost isothermal due to the absence of nuclear 
reactions. We also notice that no convective region is 
attached to the H-burning layer in these models. 

Fig. [5] shows an echelle diagram of the power spec- 
trum of HD 49385, on which we can see that the curva- 
tures of the £ = and £ — 1 ridges arc quite different. 
In fact, the distance between the £ — ridge and the 
£ = 1 rid ge is sensitive to the inner parts of the star 
(see e.g. IProvost et all |2006h ). A first attempt was 
made to find models which could reproduce the posi- 
tion of HD 49385 in the HR diagram, its large spacing, 
and the distance between the £ = and £ = 1 ridges at 
the same time. In the Main Sequence, no such model 
was found. We plotted in Fig. [6] an echelle diagram 
of the eigenfrcquencies for one of the models fitting the 
other parameters (further referred to as model A) . The 
agreement for £ = and £ = 2 ridges is satisfactory, 
but the £ = 1 ridge is very poorly reproduced. On the 
contrary, we found a PoMS model with the right T e g, 
L, Ais, and fitting the £ = 1 ridge much better (model 
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Fig. 6 Echelle diagram of the power spectrum of 
HD 49385, computed with a large spacing of Av = 56.2 /ddz. 
The frequencies of the modes are overplotted for models A 
(blue) and B (red). Squares represent I = modes, trian- 
gles, 1=1 modes and diamonds, 1 = 2 modes. A small 
frequency shift (of the order of 1 ^iHz) is applied for both 
models to partially correct for surface effects. 

B). In this model, we remark the presence of an I — 1 
avoided crossing around v = 750 /iHz, which induces a 
distortion of the ridge very similar to the one predicted 
by the analogy we presented in Sect. [5] The observed 
i = 1 ridge seems to show the same bend, which cannot 
be reproduced by MS models since they are not evolved 
enough to show t = 1 mixed modes at the appropriate 
frequency. Under the assumptions we made about the 
physics of the star, our models arc clearly in favor of a 
post Main Sequence status for HD 49385. This result 
however needs to be confirmed by further investigation, 
and a more thorough modeling of the star. 



much faster than the frequencies of the p modes. This 
is all the more true for PoMS stars, as seems to be the 
case for HD 49385, since in their interior, the increase 
of the buoyancy frequency iV is not only due to the in- 
creasing density in the center, but also to the gradient 
of chemical composition created by the burning of hy- 
drogen in shell. The observation of an avoided crossing 
in the i = 1 ridge of HD 49385 should therefore lead to 
very strong constraints on the age of the star. 

Besides, the distortion we observe in the I = 1 ridge 
is immediately related to the size of the evanescent zone 
between the cavities. This zone is limited by the turning 
point of the p-mode cavity, where u> = Si, and the turn- 
ing point of the g-mode cavity, where u> — N . The loca- 
tion of the latter turning point depends directly on the 
chemical composition in the center of the star. By ob- 
taining information on the size of the evanescent zone, 
we can therefore hope to constrain processes of trans- 
port of chemical elements in the inner parts of the star. 
For example, even though PoMS models for HD 49385 
do not have a convective core, they had one during the 
whole Main Sequence, and we expect the distortion of 
the ridge to be sensitive to the amount of overshooting 
which existed at the boundary of the core. 
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4 Conclusion and perspectives 

We presented the study of an avoided crossing involv- 
ing n modes, by extending a n analogy developed in 
Christensen-Dalsgaard ( 2003f ) to describe a two-mode 



avoided crossing. We showed that for modes of low de- 
gree i and provided the star is evolved enough, the effect 
of avoided crossings should be visible on more than two 
modes, and the presence of an avoided crossing should 
induce a characteristic and easily identifiable distortion 
of the ridge of degree I in the echelle diagram, thus 
making it easier to detect the presence of mixed modes 
in the spectrum of the star. We remark that this dis- 
tortion will be strongest for I = 1 modes, and it could 
possibly be used to identify ridges in pulsating stars 
evolved enough to show avoided crossings. 

As stated before, the frequency at which an avoided 
crossing occurs in the spectrum of an evolved star varies 
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